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1. Introduction 

An interesting aspect of nuclear dynamics is the co-existence, in atomic 
nuclei, of regular and chaotic states [1]. This is a vast subject and, for reasons 
of space, in this paper we limit ourselves to a few examples only. 

In order to highlight the difference between the regular and chaotic states, 
it is perhaps useful to remember that the low energy states (0 ~ 4 MeV 
above the ground state), termed regular, are described by a variety of models: 
the shell-model, the collective model and their various extensions [2-8]. By 
means of these models, all the properties of nuclear levels such as excitation 
energies, transition probabilities, magnetic and quadrupole momenta, etc. 
may be accurately calculated. 

When the excitation energy increases, the level density rises, making the 
calculation of nuclear properties in terms of single levels both of no physical 
interest and impossible. Instead, a statistical description, called Statistical 
Nuclear Spectroscopy (SNS), is used, based on the division into global and 
local nuclear properties [9-17]. A typical example of this is the separation of 
the level density into a global component, the secular variation, and a local 
component, the fluctuations, which are well described by the random matrix 
ensembles [18-21]. 

In recent years, the study of quantum levels in classically chaotic regions 
has shown that they have the same fluctuation properties as those predicted 
by random matrix ensembles [28-30] in a large energy range; we term these 
states chaotic. The physical significance of statistical concepts in atomic 
nuclei is therefore understood through the link with chaotic motion in hamil- 
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tonian dynamics [22-25]. 

In the first part of the present work, we review the state of the art of 
nuclear dynamics and use a schematic shell model to show how a very sim- 
ple and schematic nucleon-nucleon interaction can produce an order— ^chaos 
transition. The second part is devoted to a discussion of the wave function 
behaviour and decay of chaotic states using some simple models. 

2. Regular states in atomic nuclei: shell and collective models 

As is well known, the regular states of atomic nuclei are described, roughly 
speaking, by the mean field approximation of the shell model and by the oscil- 
lations about the mean field which give rise to collective excitations. Nuclear 
potential is essentially symmetric; nucleons move in regular orbits and, be- 
sides the energy, there are other constants of motion and other quantum 
numbers [6,7]. For the sake of completeness we report below some basic 
notions of the shell model and collective models. 

In the shell model [5] the nuclear states, like the electronic states in atoms, 
are described in terms of the motion of nucleons in a mean-field. But, as 
is well known, while the nuclear field is generated by the interactions of 
the nucleons, the atomic field is mainly governed by the interaction of the 
electrons with the nucleus. 

Experimental data suggests a shell structure for the atomic nucleus, i.e. 
the greater stability of the nuclei with 2, 8, 20, 28, 50, 82 and 126 neutrons 
or 2, 8, 20, 28, 50, 82 protons. These nuclei are called magic nuclei and their 
nucleonic numbers are termed magic numbers. The low energy levels of these 
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magic nuclei are very high, so their structure is particular stable. In the shell 
model the magic numbers represent the numbers of nucleons that saturate 
the nuclear shells. Magic nuclei are to nuclei what noble gases are to atoms. 
The hamiltonian of the shell model can be written as: 

Ho = it-J?-^ + U t , (2.1) 

where A is the number of nucleons, and Ui is the mean-field. The choice of 
the mean-field is crucial; Ui may be obtained by the usual methods using 
the many body theory. However, phenomenological approximations have 
traditionally been used instead: 

Ui = Uc,i(r) + ULS,iLi-Si, (2.2) 

where Uc,i(r) is the central potential, i.e. with a spherical symmetry, and 
Uis,iLi ■ Si is the spin-orbit interaction. It is common to use a central po- 
tential that behaves like the charge density of the nucleus (Saxon-Woods 
potential): 

= 1 + expK ^ flo)/a| . (2.3) 

where R is the nuclear radius, a ~ 0.53 • 10~ 15 m and U ~ (50-60) MeV. 
With this potential the calculations are quite complicated, so simpler poten- 
tials are usually used. 

To obtain a good agreement between the shell model results and the 
experimental data, it is necessary to add a residual interaction Hr so that 
the total hamiltonian H can be written: 

H = H + Hr, (2.4) 
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where Hr is the part of nucleon-nucleon interaction not included in Hq. 
Using second quantization formalism [8], we can write: 

A 

H = Y^ Qo+ai, (2.5a) 
i=i 

A 

H R = V ijk iafa+aia k , (2.56) 

ijkl=l 

where and V^u are the single-particle energies and the residual interac- 
tions respectively. The operators af and cij are the creation and annihilation 
operators of the ith single nucleon state: 

[a h a+]+ = 8ij, [at, aj]+ = [a+, a+] + = 0, (2.6) 

where [a, b] + = ab + ba. H and Hr can be calculated by the Hartree- 
Fock equations [8], starting from the free nucleon interactions; but often, as 
mentioned above, phenomenological approximations are used (the values of 
ti are obtained by the experimental spectrum of the nuclei with only one 
nucleon added to a double magic core). 

In order to obtain the different observables, the following procedure is 
generally used. As first step we solve the unperturbed equation: 

H \<j )a >=E^\ ( j )a >, (2.7) 

where = Y,f=i U- To obtain the eigenvalues and eigenstates of H, the 
Schrodinger equation: 

H\^ n >= E n \^ n >, (2.8) 

must be solved and \ip n > expanded by a complete unperturbed base of 
eigenstates: 

\i> n >=Y.C na \<p a > ■ (2.9) 

a 
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In this way we obtain the secular problem: 

£4 0) <W+ < <t>a\H R \<t>p > C n/3 = E n C na . (2.10) 

P 

It is clear that the sums in (2.9) and (2.10) are infinite and so the secular 
problem is not solvable. To overcome this difficulty, it is standard procedure 
to cut the basis states by introducing a finite number of configurations which 
are sufficient to describe the first excitation states. Many nucleons are frozen 
in the deeper shells of the mean field potential and form an inert core; only 
a few nucleons partially populate the single particle shells outside the core. 
These are called valence-nucleons. So there are N valence-nucleons, m active 
shells and a finite number of energy levels. 

Now, if we know the states \ip n >, it is possible to calculate all the ob- 
servables A, which characterize a nuclear state. To obtain the value of A in 
the state \ip n >, we must calculate the diagonal element < ra |A|0 ra > of the 
operator A associated to A: 

Af n =< ^n\A\^ n >=J2 C naCnf3 < <f>a\M<l>P > • (2.11) 

The non-diagonal elements < ipf\A\ipi > are associated with transition prob- 
abilities between the initial states \ipi > and the final states \ipf >■ 

If the number of valence-nucleons is high, this procedure becomes very 
complicated; on the other hand, the low energy spectrum of nuclei with many 
nucleons outside the closed shells shows a simple behaviour which changes 
systematically from one nucleus to another. These regularities are explained 
by describing the nuclear correlations with collective motions, corresponding 
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to variations in the shape of the nucleus. In this way, we obtain a general- 
ization of the shell model: the mean field is not an isotropic static potential, 
but becomes a variable field that can assume various shapes, not only with 
spherical symmetry, but also, for example, ellipsoidal. For these nuclei the 
collective motions can be divided into rotations and vibrations [6,7]. The first 
correspond to the rotation of the nuclear orientation with shape conservation 
and have low excitation energies; the second correspond to oscillations of the 
nucleus around its equilibrium shape. 

The nuclear surface in polar coordinates can be written: 

oo A 

R(9, 0) = i? {l + £ E </>)} = Ro + AR, (2.12) 

A=0/i=-A 

where Rq is the mean nuclear radius, Y\^(9, 0) are the spherical harmonics 
and axfj. are the coefficients which describe the deformation, generally time- 
dependent. Small oscillations around the equilibrium shape may be described 
by harmonic oscillators, so the hamiltonian of the system can be written: 

Hx = ^ E I«a/ + ^ E Kf, (2.13) 

fi=— A n=—X 

where the coefficients B\ and C\ are related to the frequency vibration u\ 
by: 




(2.14) 



If we consider the low energy spectrum (0 ~ 4 MeV), only A = 2, 3 are 
important, because A = represents a compression (or dilatation) without 
shape change and A = 1 represents a translation of the entire nucleus or a 
dipole oscillation. Both modes are outside the energy range considered. So 
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each quantum excitation, called phonon, has an energy E\ = hu\ and spin 

I = Xh. 

For deformed nuclei, the rotation excitation energies may be written: 

E ^ = ij^I+l)-k 2 ] + £-e, (2.15) 

where I is the total angular momentum of the nucleus and k the projection of 
I on the symmetry axis. J{= J\ = J 2 ) and J 3 are the momenta of inertia of 
the nucleus referred to the principal axes. Equation (2.15) becomes simpler 
in the case of the k=0 rotational band: 

El = 27 /(/ + 1) ' (2 - 16) 

where the allowed values for the angular momentum I are 0, 2, 4, 6, ... . 
For low angular momenta (I < 8h), there is a good agreement between the 
energy values of the rotational model and the experimental data. 

3. Chaotic states in atomic nuclei 

In the low energy excitation spectrum of an atomic nucleus, the average 
level density p(E) is small and, as discussed in the previous section, one 
might expect to be able to describe most of the states in detail using nuclear 
models. However, the average level density increases very rapidly with the 
excitation energy E (Bethe's law): 

C 



P(E) = . exp(AVE=A), (3.1) 

(E — A) 4 

where A,C,A are constants for a given nucleus [2]. Therefore, once the region 
of the neutron emission threshold is reached (~ 8 MeV), the number of levels 



is so high that a description of the individual levels has to be abandoned. 
The aim of nuclear models at this and higher excitation energies is rather 
to describe special states, like giant resonances and other collective states, 
which have a peculiar structure. But the detailed description of the sea of 
background states around the collective ones is fruitless. Thus, for example, 
observations of levels of heavy nuclei in the neutron-capture region give pre- 
cise information concerning a sequence of levels from number N to number 
(N + n), where N is an integer of the order of 10 6 (see fig. 1) [16]. For these 
densities of states a level assignment based on shell or collective models be- 
comes a very difficult task. It is therefore reasonable to inquire whether the 
highly excited states may be understood from the opposite point of view, as- 
suming as working hypothesis that the shell structure is completely washed 
out and that no quantum numbers other than spin and parity remain good. 
The outcome of such an approach is termed a statistical theory of energy 
levels [10-12]. 

The statistical theory will not predict the detailed sequence of levels in 
any one nucleus, but it will describe the general appearance and the degree 
of irregularity of the level structure that is expected to occur in the atomic 
nuclei, which are too complicated to be understood in detail. 

We describe a complex nucleus as a "black box", in which a large number 
of particles are interacting according to unknown laws [13]. The problem is 
then to define, in a mathematically precise way, an ensemble in which all 
possible laws of interaction are equally probable. This program, initiated 
by Wigner and developed by many authors, has, to a large extent, been 
successful [13] . 
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3.1 The random matrix theory 



An appropriate means to define an ensemble of random matrices is pro- 
vided by the random matrix theory [14,15]. The hamiltonian matrix H is an 
NxN stochastic matrix (its matrix elements hij are random variables) and 
the probability density is specified by: 

P(H)dH = P(h u , h 12 , h NN )dhudh 12 ...dh NN , (3.2) 

where the probability information is: 

I{P(H)} = J dHP(H)\nP N (H). (3.3) 

The aim is to find the function P(H) that minimizes I. This is equivalent 
to assuming the least possible knowledge about P{H). We impose that the 
are real and, to limit the eigenvalues of H to a finite range, a condition 
on its norm [Tr(if 2 )] 1//2 is also imposed. Thence P(H) should minimize I 
subject to the constraints: 

J dHP(H) = l,j dHP(H)Tr(H 2 ) = C, (3.4) 

which leads to: 

P(H) = exp {Ax + \ 2 Tr(H 2 )}. (3.5) 

By inserting (3.5) in (3.4), the Lagrange multipliers Ai,A2 may be deter- 
mined. This probability distribution defines the so-called Gaussian Orthog- 
onal Ensemble (GOE), since P(H) is invariant for orthogonal transformation 
and the elements h^ are independent random variables. We can observe that, 
without the time invariance, the are complex numbers and we have the 
Gaussian Unitary Ensemble [14]. 

10 



With the help of (3.5) one can easily obtain P(E) = P(E\, E%, En) 
and so the statistic of the energy spectrum. Mehta defines "statistic" thus: 
U A statistic is a quantity which can be calculated from an observed sequence of 
levels alone, without other information and whose average value and variance 
are known from the theoretical model. A suitable statistic is one which is 
sensitive for the property to be compared or distinguished and is insensitive 
for other details" [15]. Various statistics may be used to show the local 
correlations of the energy levels; in this work we shall use P(s) and A 3 (L) 
mainly. 

The first statistic measures the probability that two neighboring eigen- 
values are a distance "s" apart, in the average level distance unit. For this 
statistic the GOE average is closely approximated by the Wigner distribu- 
tion: 



which gives level repulsion. 

The second statistic is defined for a fixed interval (—L/2, L/2), as the 
least-square deviation of the staircase function N(E) from the best straight 
line fitting it: 



where N(E) is the number of levels between E and zero for positive energy, 
between — E and zero for negative energy. The A 3 (L) statistic provides a 
measure of the degree of rigidity of the spectrum: for a given interval L, 
the smaller A 3 (L) is, the stronger is the rigidity, signifying the long-range 



P(s) = -sexp(--s ), 



(3.6) 




L/2 



(3.7) 
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correlations between levels. For this statistic in the GOE ensemble: 



ML) 



(3.8) 



. ^logL, L>1 

For the GUE there are similar results [14]. 

If the mean level density for the GOE is calculated, we obtain: 



( 



p(E) 



V2Na 2 - E, \E\ < 2Na 2 



(3.9) 



. 0, \E\ < 2Na 2 

Although this is an unrealistic result, it can be explained by remembering 
that global and local behaviours are on different scales and GOE is a good 
model for local properties. 

3.2 Comparison of GOE predictions with experimental data 

Neutron resonance spectroscopy on a heavy even-even nucleus typically 
leads to the identification of about 150 to 170 s-wave resonances with J n = 
| + located 8-10 MeV above the ground state of the compound system, with 
average spacings around 10 eV and average total widths around 1 eV. Proton 
resonance spectroscopy yields somewhat shorter sequences of levels with fixed 
spin and parity, with typically 60 to 80 members. For the statistical analysis, 
it is essential that the sequences be pure (no admixture of levels with different 
spin or parity) and complete (no missing levels) [17]. Only such sequences 
were considered by Haq, Pandey and Bohigas [18]. Scaling each sequence to 
the same average level spacing and lumping together all sequences one leads 
to the "Nuclear Data Ensemble" (NDE), which contains 1726 level spacings. 

12 



As shown in fig. 2 the agreement between the experimental data and the 
GOE predictions is surprisingly good (in the GOE model there are no free 
parameters) [19]. 

It is also important to observe that this agreement in not limited to 
nuclear levels, but is valid for all quantum systems (nuclear, atomic and 
molecular) in a high energy range [1,29]. 

The behaviour of spectral statistics near the ground state is also of con- 
siderable interest. Garrett, German, Courtney and Espino [20] (fig. 3) and 
also Shriner, Mitchell and Von Egidy [21] (fig. 4) have shown that for these 
energies nuclei do not follow the GOE results but behave like a Poisson en- 
semble or with an intermediate behaviour between Poisson and GOE. In the 
Poisson ensemble: 



which gives a high probability for the occurrence of near degeneracies, no 
level repulsion and no correlation between levels. 

To sum up, we have seen that the spectra of atomic nuclei show a Poisson^GOE 
transition on increasing the excitation energy (see fig. 5) over the yrast line. 
In the next section we show how this behaviour is not of a purely quantal 
nature but has a classical counterpart in the regular and chaotic hamiltonian 
systems. 

4. Order and chaos in classical and quantum mechanics 

In this section we introduce the dynamical systems and their stability, 
which may be studied by means of the Lyapunov exponents and metric en- 





L 

15' 



(3.10) 
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tropy. With the help of these quantities, we clarify the concept of ergodic 
system giving a hierarchy of chaos. Then we extend the study to hamiltonian 
systems, which are essential in order to study the transition order— ^chaos in 
nuclear physics. 

A dynamical system [22-25] is defined by N differential equations of the 
first order: 

j t z(t) = f(m,t), (4.i) 

where the variables z— (z±, ...,zn) are in the phase space Q (the euclidean 
space R N , unless otherwise specified). These equations describe the time 
evolution of the variables and the system they represent. 

A solution of the dynamical system is a vector function z(z ,t), that 
satisfies (4.1) and the initial condition: 

z(z , 0) = z , (4.2) 

often written simply z(t) without the initial condition dependence. 

The time evolution of z G Q is obtained with the one parameter group of 
diffeomorfism g l : Q — > Q, so that: 

±(gtz)\ t=0 = f(z,0). (4.3) 

The group g l is called phase flux and the solution is called orbit. The system 
is called hamiltonian, if the dimension of Q is even and there exists a function 
H(z, t) given by: 

f(z(t),t) = JVH(z,t), (4.4) 

where: 

(4 - 5 » 
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is the symplectic matrix. H(z,t) is the hamiltonian. 

On the phase space Q one usually defines a probability measure fj, : Q — > 
f2, so that = 1. If we choose a subspace A of f2, the system is measure 
preserving if: 

ptfA) = p(A). (4.6) 

It is well known that hamiltonian systems preserve their measure: the Liou- 
ville measure. Dynamical systems which do not preserve their measure are 
called dissipative, and usually have a measure contraction in time evolution. 

The dynamic of a system is called regular if the orbits are stable to in- 
finitesimal variations of initial conditions. It is called chaotic if the orbits 
are unstable to infinitesimal variations of initial conditions [23-25]. Useful 
quantities to calculate this behaviour are the Lyapunov exponents, which 
give the stability of a single orbit, and the metric entropy, which represents 
a mean exponent for the entire system. A vector of the tangent space TQg 
to the phase space fl in the position z is given by: 

u{z)=\ ^ S) -^\ (4.7) 

where q(0) = z and q(s) G fl The tangent space vectors are the velocity 
vectors of the curves on M; there are obviously N independent vectors. 
Now we can define the Lyapunov exponent: 

X(z) = lim - \n\cj(t)\, (4.8) 

t^oo t 

where u(t) is a tangent vector to z(t) with the condition that |tJ(0)| = 1. 

It can be demonstrated that the limit given by (4.8) exists for a compact 
phase space, and that it is metric independent. Fixing an orbit in the N 
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dimensional phase space, there are N distinct exponents Ai, A at, called first 
order Lyapunov exponents. If the orbit has positive Lyapunov exponents, it 
is chaotic. 

To characterize globally the chaoticity of a system, we introduce the met- 
ric entropy. If a(0) is a partition of non-overlapping sets that completely 
cover the phase space Q at the initial instant: 

a(0) = {A(0) : E(E) = U^(0), A(0) n ^(0) = 0}, (4.9) 

the partition a(0) can evolve in a discretized time flux: 

a(0),a(l),a(2),...,a(n). (4.10) 

Then we define f3(n) as the intersection set of all the sets at every instant: 

(5{n) = {B,(ra) : £,(n) = A(0) n A,-(l) n ... n A k (n)}. (4.11) 

In this way, the number of sets defined by f3(n) do not decrease when n is 
increased. By introducing a probability measure \i preserved by dynamics, 
one can define the metric information of the partition Bin): 

^ = -Y, l i(B i (n))]n»(B i (n)). (4.12) 

i 

The information I n has the following properties: 

(i) I n — if, and only if, there exists a set B^n) of the partition f3(n) so 
that fx(Bi(n)) = 1; 

(ii) I n assumes the maximum value when /3(n) = {Bi(n), ...,B^(n)} and 
in this way fj,(Bi(n)) = jj Vi i.e. /„ = IniV. 
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The metric entropy of the partition a(0) is defined as: 

h(AM,fi) = lim — (4.13) 

i.e. a time average of the metric information. If all the elements of a(t) have 
a measure that decreases, on average, exponentially with t, then the metric 
entropy will be positive. A positive entropy indicates that there does not 
exist a finite number of measures to guess the next one. 

We indicate the maximum value of h on all partitions Ai with hxs' 

h KS = max{/i(A(0),//),VA(0)}, (4.14) 

called Kolmogorov-Sinai entropy. According to a very important theorem 
[26]: 

h KS (n) = [ dfi(z) ]T \ t (z) (4.15) 

with A subspace of Q and Aj Lyapunov exponents. Therefore the Kolmogorov- 
Sinai entropy is a very useful tool for showing chaotic behaviour in the region 
A. 

A system is called ergodic if the time average is equal to phase space 
average: 

lim - f dtf^zit)) = I d»(z)f(z). (4.16) 
t~*oo t Jo Jn 

Incidentally, as is well known, Boltzmann started from the "ergodic hypoth- 
esis" to obtain statistical mechanics of equilibrium. But ergodicity is not 
sufficient to reach an equilibrium state: one must consider mixing systems. 

In a mixing system, every finite element of the phase space occupies for 
t — > oo the entire phase space f2; more precisely: V7L, BcH with fi(A) and 
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KB) ^ 0, 

n(B H g t A) , AS .. 
fe - ,m =M^- (4-17) 

To have quantitative information of orbit separations, we must introduce 

K-systems (Kolmogorov), which are mixing systems with a positive metric 

entropy: 

h KS > 0. (4.19) 
Such systems are typical chaotic systems. 

Among the K-systems, the most unpredictable ones are the B-systems 
(Bernoulli), which have the Kolmogorov-Sinai entropy equal to the entropy 
of every partition: 

h KS = h(A l (0), f i), VA(0). (4.20) 
4.1 Classical chaos 

In classical mechanics, the state of a system of coordinates and moments 
Pi, i = l,...,n in the N=2n dimensional phase space Q, is specified by the 
hamiltonian H = H(p,q), with p = (pi, ...,p„), q= (q u ...,q n ) [22]. 

As is well known, the time evolution is obtained by the Hamilton equa- 
tions: 

9H . OH , t . 

Qi=~-, Pi = -o~- 4 - 21 

dpi dqi 

These equations, with the position z— (q 1 , ...,q n ,pi, ...,p n ), can be written 
in the more compact form (4.4). 

The hamiltonian system is integrable if there are n functions defined on 

SI: 

F l = F t (z) i — 1, N (4.22) 
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in involution: 

n BF BF BF BF 

l F *> F Ar* = T,PP--p-7T = <>> V ^ ( 4 - 23 ) 
t^i oq k op k dq k dp k 

and linearly independent. [ , ]pb are the Poisson Brackets. 
For the conservative systems we have F\ = H{z) and also: 

^ = [H,F t ] PB = 0. (4.24) 

Because there are n constants of motion, every orbit can explore only the n 
dimensional manifold Qf m . 

Q f = {z:F l (z)=f t ,t=l,...,N} (4.25) 

If flf is compact and connected, it is equivalent to a n dimensional torus: 

T n = {(0i, ...,0 n ) mod 2tt} (4.26) 

There are n irreducible and independent circuits 7$ on Qf and there exists a 
canonical transformation: 

$ (4.26) 
generated by the function S(q,I), so that: 

I i = j>dq-p 1 Oi=Qj- (4-27) 
The ij are called action variables and the 6i are called angle variables. 

— * 

The moments p and coordinates q are periodic functions of 9 with period 2tt: 
5 = E ^('V™'*, P = £Pm('V™'* ( 4 -28) 
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where m = (mi, ...,m n ) is an integer vector. 

The hamiltonian depends only on action variables: H = H(l), and so the 
new equation of motion are: 

so: 

0i = ^(7)* + ^ (4.30) 

where a, are constants and cjj are the frequencies on the torus. 

The orbits of an integrable system are quasi-periodic, with n quasi- 
periods: 

Ti = — , (4.31) 

UJi 

and the orbit is closed if there exists a period r so that: 

0(r) = (9(0) + 2tt/3, (4.32) 

— * 

with D integer vector. To have the closure, we must have: 

u ■ rhi = 0, (4.33) 

and the period is: 

2tt 2ttA A 

r = — = = — -. (4.34) 

LO c UJi Ti 

If we do not have n independent relations among frequencies, the orbits do 
non close: we have the so called "torus ergodicity", but the motion is still 
regular (see fig. 6). 

Adding a small perturbation V(I,9) to an integrable hamiltonian H (I), 
the total hamiltonian can be written: 

H(I, 6) = H (l) + X V(I, 9), (4.35) 
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and, generically, the integrability is destroyed. As a consequence, parts of 
phase space become filled with chaotic orbits, while in other parts the toroidal 
surfaces of the integrable system are deformed but not destroyed; thus we 
have a quasi-integrable system. 

With growing x, chaotic motion develops near the regions of phase space 
where all the u>i are commensurate; or, more precisely, where the equation 
(4.33) is obeyed by integer vectors. 

Conversely, tori of the integrable system on which the Ui are incommen- 
surate are deformed, but not destroyed immediately (KAM theorem). As x 
increases, the phase space generically develops a highly complex structure, 
with islands of regular motion (filled with quasi-periodic orbits) interspersed 
in regions of chaotic motion, but containing in turn more regions of chaos. As 
X grows further, the fraction of phase space filled with chaotic orbits grows 
until it reaches unity as the last KAM surface is destroyed. Then the motion 
is completely chaotic everywhere, except possibly for isolated periodic orbits 
[22,23]. 

It is very useful to plot a 2n— 1 surface of section Pcfl, called Poincare 
section (see fig. 7). As shown in figure 8, for an integrable system with 
two degrees of freedom, the x = Poincare section of a rational (resonant) 
torus is a finite number of points along a closed curve, while the section of 
an irrational (non resonant) torus is a continuous closed curve. 

Adding a perturbation, the section presents closed curves (KAM tori), 
whose points are stable (elliptic), and also curves formed by substructures, 
residua of resonant tori, whose points are unstable (hyperbolic). As the 
perturbation parameter increases, the closed curves are distorted and reduced 
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in number. 



4.2 Quantum chaos 

In non relativistic quantum mechanics, the time evolution of a state \ip > 
of a system with hamiltonian H is given by the Schrodinger equation [27]: 

ih—\^>=H\^j>. (4.36) 

The unitary time-evolution operator U (t) is: 

tty)=exp(~), (4.37) 

so that \ip(t) >= U(t)\ip >. If we choose a new state j^' > near to \ip > we 
have: 

< ^'(OM*) >=< ^'|f> + (t)f>(t)|^(t) >=< V>V > • (4.38) 

The distance between two quantum states is always time- independent, so 
there is no chaos in quantum mechanics in the sense of time exponential 
divergence of near states. However an interesting question is: whether there 
are quantum-mechanical manifestations of classical chaotic motion. We shall 
use the term quantum chaotic system in the precise, and restricted, sense of 
a quantum system whose classical analogue is chaotic [24,25]. 

The studies of the eigenvalues of billiards with different shapes [28,29] 
have shown that if the system is classically integrable the spectral statistics 
are well modelled by the Poisson ensemble (3.10) and if the system is clas- 
sically chaotic, by the GOE (3.6, 3.8) or GUE ensembles, depending on the 
time-reversal symmetry (fig. 9). 
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When the classical dynamics of a physical system is regular, the short- 
range properties of the corresponding quantal spectrum tend to resemble 
those of a spectrum of randomly distributed numbers (the Poisson spec- 
trum). This is because regular classical motion is associated with integrabil- 
ity or separability of the classical equations of motion. In quantum mechanics 
the separability corresponds to a number of independent conserved quanti- 
ties (such as angular momentum), and each energy level can be characterised 
by the associated quantum numbers. Superimposing the terms arising from 
the various quantum numbers, a spectrum is generated like that of random 
numbers, at least over short intervals. When the classical dynamics of a phys- 
ical system is chaotic, the system cannot be integrable and there must be 
fewer constants of motion than degrees of freedom. Quantum mechanically 
this means that once all good quantum numbers due to obvious symmetries 
etc. are accounted for, the energy levels cannot simply be labelled by quan- 
tum numbers associated with certain constants of motion. The short-range 
properties of the energy spectrum then tend to resemble those of eigenvalue 
spectra of matrices with randomly chosen elements. 

Berry and Tabor [30] predict level clustering for any integrable system 
(with the exception of the harmonic oscillator that has equal spaced levels 
[31]) in the asymptotic high energy regime, using the semiclassical {Ti — > 0) 
Einstein-Brillowin-Keller (EKB) quantization of action variables [24]: 

E ni ,...,n N = H(h = him + ai/4), ...,I N = Hn N + ajv/4)), (4.39) 

where a; is the Maslov index for the J; action variable. In one dimension: 
a = for the rotator and a — \ for the oscillator. 
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In quasi-integrable and chaotic systems the EKB quantization is not ap- 
plicable but we can calculate the level density p(E) using a semiclassical 
formula obtained by Gutzwiller [32,33]: 



p(E)=p(E)+ Posc (E), (4.40) 



where: 



?{E) = J^Kf / im(E ~ m (441) 
p „ c(B) ^_2_ cos( SM_ ?) . (4 , 2) 

The summation is over all the periodic orbits of the classical phase space, Tj 
is the j-periodic orbit period, / is a function of the j-periodic orbits Lyapunov 
exponent and Sj(E) = §^ dq- p. 

With the help of (4.42) Berry [34] has calculated the semiclassical be- 
haviour of the spectral rigidity A^(L) for a chaotic time-reversal system: 



A 3 (L) 



4logL, URU , (4.43) 



^ 2^2 \ogeL max + oscillations, L > L max 
where L max = with T min the period of the shortest periodic orbits of 

J- rain 

the classical phase space. The long-range periodic orbits contribute to the 
universal behaviour, but there is also a non universal behaviour not predicted 
by GOE. 

In the next section, we describe a schematic shell model to show how an 
order to chaos transition may occur by studying the classical and quantum 
properties of the model. 
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5. Order to chaos transition in a schematic shell model 



To study the transition from regular to chaotic states in atomic nuclei, 
an extension of the two level Lipkin-Meshkov-Glick (LMG) model was pro- 
posed by Meredith, Koonin and Zirnbauer (MKZ) [35-37]. The two level 
LMG model has only one degree of freedom, i.e. the particle number in the 
upper level, and so its classical limit does not have chaotic behaviour. The 
generalization proposed by MKZ has the SU(3) symmetry and two degrees 
of freedom: its classical limit shows an order to chaos transition. 

The SU(3) model consists of M identical particles, labelled by the index 
n, each of which can be in three single-particle states having energy e^. The 
hamiltonian is: 

£ = E*£«+E^-g?,-> (5-1) 

i=l ij=l 

with e 3 = e x = e and e 2 = 0, where is the strength of the two-body 
interaction between states i and j: 

Vi^Vji and Vu = 0, (5.2) 

with: 

M 

Gij = E V-nfinj- (5-3) 
n=l 

The operator d^ creates a particle n in state i, and d n i annihilates a particle 
n from state i; each term of takes a particle, n, out of state j and into 
state i, and each Ga counts the particles in state i. 

Applying the usual fermion anti-commutation rules to the d ni operators, 
we get: 

[G^, Ghk]- = GikSjh — GjhSik, (5.4) 
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and the 9 operators Gij are generators for U(3). Taking into account the 
number conservation, Ya=i Gu = M, this becomes SU(3). 

5.1 Classical limit 

The quantum 4-dimensional phase space of SU(3) is U(3)/U(2) <g> U(l), 
where C/(2)<g)[/(l) is the maximal stability subgroup which leaves the ground 
state 1 00 > invariant up to a phase factor [38]. 

The coherent states of U(3)/U(2) ® U(l) are: 

\SU(3),$ >= $|00 >, (5.5) 

where 

^ = exp{X)(i/ i G i i-^G li )}, (5.6) 

i=2 

with Vi complex numbers. 
The Bergmann kernel is: 

k(z, z*) =< SU{3), &\SU(3), $ >= (1 + Z + Z) M , (5.7) 

where Z is the 3x1 matrix: 

Z = fo) = M *S&1, (5.8) 

with v = y/iwl + v z vl. The metric of U(3)/U(2) <g> U(l) is then given by: 

8 2 lnfc(z,z*) 

= dzM ' (5 ' 9) 

and the nondegenerate closed two-form u of U(3)/U(2) <g> U(l) is: 

3 
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with Poisson Brackets: 



[/ '^ PB = -\f/ (5 - U) 
By introducing the canonical coordinates: 

-4= (ft + = — sin z/, (5.12) 

one can show that (5.12) will be transformed into the following canonical 
form: 

v>fo°=h%-g-g-%-\- ("3) 

■j±2 d 1i 9 Pj d Pj d( li 

The semiquantal dynamics of this many-body system with fixed nucleons 
number M is determined by varying the effective action: 

S = J [\{p ■ dq- q- dp) - Hip, q)dt}. (5.14) 

The result of such a variation is the set of classical-like dynamical equations: 

dqj _ dHjp, q) dp t _ dHjp, q) 

dt dpi ' dt dqi 

with the Hamiltonian functions given by: 



(5.15) 



H{p, q) =< SU(3), $\H\SU(3), $ > . (5.16) 

Equation (5.15) is in fact equivalent to the time-dependent Hartree-Fock 
(TDHF) dynamical equations. It should be noted that (5.15) is not a classical 
limit because there is at least first-order quantum correlation in Hip, q). 

The phase space representation of the generators in the coherent- 
states basis is: 

< |G n | >= l -i2M-p 2 -q 2 ), (5.17) 
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<\Gn\ >= g (ft + ^) V / 2M- P 2_ g 2 ! 

< >=< >*, 

where | > represents \SU(3), $ > and (g^ + ?! + pl + pD < 2M. 

The quantum correlations of the quadratic function of the generators can 
be computed and the results are: 

A(GJ-) =< |<%| > - < |G y | > 2 = < 14-1 > 2 • (5-18) 

In this way the mean value of the hamiltonian in the coherent states repre- 
sentation is: 

H(P, 0) = |[2 - (^ +Pl + ql + ql)} + f (p£ + <7 2 2 ) + |(p§ + <7 3 2 ) + (5-19) 

+ f [1 " ^] [(P2 + PS)' - («2 + ft 2 ) 2 + (?2 - PlM ~ Pl) + %2<73P2P3)]. 

It should be noted that the phase space has been scaled in (5.19) such that 
ql + ql + p\ + P3 < 2, with x = MV/e. The classical limit of the SU(3) 
model can be realized by the limit M — > 00, and, in this limit, the quantum 
correlation A(G?) goes to zero (see 5.18). So the classical hamiltonian is 
(6=1): 

Hci = -1 + \ql{l ~ x) + ^3(2 - X) + ^2(1 + X) + ^P 2 (2 + x) + 

+^x[(? 2 2 + <? 3 2 ) 2 - {pl + pi) 2 - {ql - pl){ql - pi) - ^qmpzl (5.20) 

with the phase space given by Q = {(?2, 93, P2, P3) £ -R 4 : {ql + ql + p\ + 
P 2 ) < 2}. This hamiltonian represents two oscillators coupled non-linearly 
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by the parameter x (the interaction between the nucleons). The system is 
integrable if the interactions are neglected but the interactions break this 
regular behaviour and produce an order to chaos transition. 



5.2 Classical calculations 

Using the Hamilton equations of (5.20), in order to analyze the stability 
of the system, we calculated the periodic orbits of this model: 

'z = JVH d (z, X ), (5.21) 

where z = (q 2 , gs,P2,Ps) , V = (^-, ^), e J is the 4x4 symplectic 

matrix: 




(5.22) 



where I is the 2x2 identity matrix. 

To explore the phase space f2, we chose a 4-dimensional lattice of initial 
conditions z = z(t = 0), z e fl The side of the lattice has about 10 4 initial 
conditions. 

The time evolution of (5.21), with the initial condition z , is obtained 
using a A th order Runge-Kutta method [39]. The following function can then 
be constructed: 

d(z ,t) = d(z ,z t ) = \z t - z \, (5.23) 

where z t = z(t). (5.23) is the function to minimize, with minimum value 
equal to 0. To obtain a periodic orbit, it is not actually necessary to vary all 
five parameters (z ,t)] in fact we can fix E = H(z ) or t, which correspond 
to focusing either on constant energy or constant period in the E-T plot. 
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Another algorithm [40] has been used to find the minimum of d. Since d 
is the result of a complicated calculation and its derivatives are not available, 
this algorithm is particularly useful in our case because the derivatives of d 
are not required [41]. 

These periodic trajectories occur in one parameter families, each of which 
describes a continuous curve in the energy-period plane (fig. 10); just as 
in quantum mechanics the set of energy levels is characteristic of a given 
hamiltonian H, so H d can be characterized by the plots E-T, where E is the 
energy of a periodic trajectory and T the corresponding period. 

If w is a vector of the tangent space TVlg of the phase space manifold Q 
at z, its time evolution is given by: 

(5.24) 

By (5.21) and (5.24) the Lyapunov exponents can be calculated: 

\{z) = Urn -In \w{t)\. (5.25) 
In terms of the stability matrix M(0, t) [42], defined as: 

*WM) = gg, (5.26) 

\{z) can be written: 

X(z) = lim -ln|M(0,0l, (5.27) 

t— >oo f 

where |M(0, t)\ is a norm of the matrix M(0, t). 

This matrix can be calculated by solving its equations of motion: 

d 2 H(z) 

M = J—±±M, (5.28a) 
oz 2 
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with the initial conditions: 

M (0) = / (5.286) 

where I is the 4x4 identity matrix. The calculation of the Lyapunov expo- 
nents is related to that of the eigenvalues U{ of the matrix M(0,T): 

\ i (z) = ^lna i . (5.29) 

Now, using the unitariety of M, a periodic orbit is unstable if 

Tr(M) > 4 or Tr(M) < (5.30a) 

and stable if 

< Tr(M) < 4. (5.316) 

It is interesting to study the change of stability of periodic trajectories as 
a function of x- 111 figure 11 the ratio between the number of stable orbits 
and the number of total orbits with period T < 30 is plotted vs x- As can be 
seen the sensitivity of the orbits to a small change of x is quite different for 
different values of x, reflecting the transition order — ■> chaos as the coupling 
constant increases. 

As mentioned in section 4.2, the calculation of the long-range behaviour 
(non-universal) of A 3 requires the knowledge of the shortest period T min of 
the closed trajectories. For a chaotic system with time-reversal symmetry 
[34]: 



— In (eL max ) — — 

where: 



A™ = — \n(eL max )--, (5.32a) 



Lmax — rr ^ (5.326) 



T 

- 1 mm 
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and p is the average level density. Obviously p and T min are functions of the 
coupling constant X- 

Figure 12 shows A 3 (L) for x = 100, T min = 0.04, p = 115. As can be 
seen, the breaking of universality is near L=38. Obviously the semi-classical 
estimate of the A 3 (L) saturation does not agree with the quantal calculations 
for x = 0.5 and x = 2 because the equation (5.32) is valid only for a chaotic 
system [43]. 

Incidentally, Aberg [44] gives a rough estimate of the shortest periodic 
orbit in a nucleus: 

T m in — — J—^ 5 (5.33) 

vf 0.3c 

and thus L max ~ 80A-^p. For l52 Dy at E x ~ 3.0 MeV, p ~ 560 MeV" 1 , 
we get L max ~ 8000, i.e. we would need a sequence of more than 8000 
energy levels with the same spin and parity. This example shows that we 
should hardly expect to see a non-generic behaviour of long-range spectrum 
fluctuations due to short periodic orbits in nuclei. 

5.3 Quantum calculations 

In order to find the eigenstates and eigenvalues of the hamiltonian (5.1), 
we need to find a complete set of basis states. A natural basis can be written 
\bc > meaning b particles in the second level, c in the third and, of course, 
M — b — c in the first level; in this way [00 > is the ground state with all the 
particles in the lowest level [36,37]. We can write the general basis state: 



\bc>= X^G 6 21 Gy00>, (5.34) 
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with y the normalizing constant. 

From the commutation relation (5.4) we can calculate the expectation 
values of || and thus, eigenvalues and eigenstates of j|; in this way the 
energy spectrum range is independent of the number of the particles: 

where: 



Q b >c>, bc = y/b(b-l)(M-b- 


-c+l)(M-b- 


- c + 2)5 b _ 2 b >5 cc > 


(5.36) 


+y/(b+l)(b + 2)(M- 


-b-c)(M-b- 


-c-l)5 b+2b ,5 ccl 




+y/c{c-l)(M-b- 


c+l)(M-b- 


c + 2 )V<W' 





+ y/(c+l)(c + 2)(M-b-c)(M-b-c- 1) V 5 c + 2, c > 



+ ^ + l)(& + 2)c(c-l)^ 5 C _^ C , 

+ v /6(6-l)( c+ l)( c + 2K_ 2 ^ 5 c+2iC , 

and x — MV/e. The expectation value < ^ > is real and symmetric. For 
any given number of particles M, we can set up the complete basis state, 
write down the matrix elements of < || > and then diagonalize < || > 
to find its eigenvalues. < || > connects only states with A6 = —2, 0, 2 
and Ac = —2, 0, 2 which makes the problem easier. We group states with 
b,c even; b,c odd; b even and c odd; b odd and c even. This means that 
< jj > becomes block diagonal containing 4 blocks which can be diagonalized 
separately; these matrices are referred to as ee, oo, oe and eo. 

To separate regular and chaotic state in a quantum system a powerful 
tool is the study of spectral statistics. We obtain [45] a good agreement with 
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GOE in the classically chaotic region and, in the classically regular region, 
a good agreement with Poisson statistic (see fig. 13). In this figure the 
continuous curve is the Brody distribution, discussed in section 6.5. 

On the basis of the semiclassical torus quantization, the presence of cross- 
ings in a small x neighbourhood may be interpreted as the signature of quasi 
crossing in the true system, and thus as a signature of torus destruction, 
when the exact levels are "split" at the crossing [46-48]. 

Another method to study the quantum stochasticity is the sensibility of 
the energy levels to variations of the perturbation parameter; the behaviour 
of the curvature of the energy level E(x) in a small x neighbourhood is an 
example [49]. 

When the parameter x — 0, then the hamiltonian represents two oscilla- 
tors, and there are many degenerations, but for x these degenerations 
are broken. We calculated [45], for a large number of particles (semiclassi- 
cal limit) the density of quasi crossings outside the degeneration region as 
function of the parameter x- 

P(X) = (5.37) 

where AiV is the number of quasi crossings in the parameter range A%. To 
obtain AiV, we fixed three values x ~ ^-Xi X an d X + ^X an d imposed that: 

Si ( X - A X ) > Si( X ) (5.38a) 

Sl ( X + A X )> Si (x) (5.386) 

where Sj(x) = E i+ i(x) — -^i(x)- The results (fig. 14) show for all the classes 
a maximum of quasi crossings for x = 2 in agreement with the transition to 
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chaos of fig. 12. 

In order to study the sensitivity of the energy levels to small changes of 
the parameter x, we have used the statistic A 2 (E) [50], defined as: 

A 2 (^) = \Ei( X + A X ) + Ei( X - Ax) - 2E i ( X )\. (5.39) 

This statistic measures the curvature of E i in a small range A%. Fig. 15 
shows A 2 (E) in the [-0.5,0.5] energy range for different values of x; in this 
case we note the formation of a peak in correspondence to the x = 2 region. 
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6. Wave function behaviour and EM decay of regular and 
chaotic states 



As we have seen in the previous sections, it is now well known that in 
atomic nuclei there are two different states: regular and chaotic. These 
states can be highlighted with the aid of various statistics. The statistics 
we have considered before concern the eigenvalues of the system, and the 
most frequently used are those of level spacing; i.e. the nearest neighbour 
distribution, P(s), and the stiffness of the spectrum A 3 (L). Since these 
statistics display clearly different behaviour depending on the regularity or 
chaoticity of the nucleus, they are very useful to distinguish between the two 
regimes and, at the same time, characterize the transition from one region 
to the other. 

However, it is well known that, in general, the wave functions and the 
transition probabilities are much more sensitive to the purity of states than 
the eigenvalues are. Consequently, it is interesting to undertake the study of 
statistics directly related to the wave functions with the aim of characterizing 
the ordered or chaotic behaviour of atomic nuclei. 

Using the SU(3) schematic shell model for the nucleus, described in Sec- 
tion 5, we have studied the various statistics concerning the wave functions, 
such as the intensity of the momenta of the wave function I m , the information 
measure or entropy S, and the correlation functions K n>rn . 

The momenta I m and the information measure S are defined [51]: 





b,c 



36 



S = -]T |^(6, c) | 2 In |^(6, c) | 2 . (6.2) 

Figure 16 shows that the lower momenta (m=2,3,4) diverge in the regular 
region but have a constant behaviour in the chaotic region, as predicted by 
[52]. 

The information measure S also shows different behaviours in the two 
regions, assuming high values in the chaotic region and low values in the 
regular region because wave functions are localized (see fig. 17). 

Another way to separate regular from chaotic states is with the help 
of the transition probabilities among different states [53]. Only a few strong 
transitions between regular states and many weak transitions between chaotic 
states are expected. The transition probability is defined: 

Pij ] = I < Vi\T k \Vj > I 2 , (6-3) 

where: 

T 1 = ^(G 12 -G 21 ), T 2 = i(G 13 -G 3 i), T 3 = i(G 23 - G 32 ). (6.4) 

|#>=XXJ&,c> (6.5) 

a=l 

is the wave function of the j-state of class r with dimension d T (ee, oo, eo 

pk 

and oe). Scaling Pf, with a local mean value < Pjj >, the distributions < pl > 
follow the above previsions very well (see fig. 18). 

The previous results clearly show that the wave function behaviour of 
atomic nuclei is a very powerful tool for distinguishing regular and chaotic 
regimes. In particular, the study of the transition probabilities between dif- 
ferent states is of great interest. 
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6.1 Decay of regular and chaotic states 

Recently, new interest has been shown in the study of the decay of unsta- 
ble quantum states. In this section we briefly review the literature of major 
interest. The section is organized as follows: first, we look at the current 
experimental stage to illustrate the reasons and interest for studying the de- 
cay of quantum systems. Secondly, we review the main current theoretical 
approaches to the problem. We begin with the study of "warm" rotational 
nuclei, to continue with that of unstable systems described by means of an 
anti-hermitian effective hamiltonian. Finally, we briefly discuss the Interact- 
ing Boson Model (IBM) that affords a wide spectrum of applications, from 
vibrational to rotational nuclei as well as transitional regions. 

6.2 Experimental stage 

The new generation of experimental instruments like GASP, EUROBALL, 
EUROGAM, etc [54], which provide a resolution in the order of a hundred 
times greater than that of the previous generation, will open a new way in 
nuclear spectroscopy research. Over the last few years, as discussed above, 
great interest has been shown in the study of the properties of the atomic 
nuclei from a statistical point of view. This has made it possible, with the 
aid of the Random Matrix Theory [14], to distinguish between regular and 
chaotic regimes, the onset of chaos, the transition from one regime to the 
other, etc. 

Let us now turn our attention to the recent statistical analysis of 7-ray 
spectra of "warm" rotational nuclei by Garrett, Hagemann and Herskind 

38 



[55]. In their work, two and three- fold energy correlation measurements 
between gamma rays were analyzed to investigate high spin states of the so 
called quasi-continuum that spreads up to 5-10 MeV above the yrast line. 
This study had already been proposed by Guhr and Weidenmuller [56] who 
stressed the coexistence of collectivity and chaos in nuclei at several MeV of 
excitation energy; they suggested the study of stretched collective E2 decays 
of high spin states in deformed even-even nuclei at, or above, the yrast line. 
It was shown that the collective E2 transitions from states of spin J, located 
several MeV above the yrast line, populate many states (of the same parity) 
of (J — 2) spin rather than a single one (as happens with the E2 transitions 
from states near the yrast line, where a well defined rotational band structure 
exists). Thus the transition strength is spread over an energy interval that 
is called rotational damping width. The analysis of the group of Herskind is 
indeed oriented in this direction. By performing a fluctuation analysis of the 
2-D E ll yiE l2 spectrum they obtained information on the number of decay 
routes the nuclear decay flow takes to go from the initial high spin state to 
the damped region (which is supposed to be in the diagonal valley, x — y, of 
the 2-D spectra). A schematic illustration of the average flow of 7-ray decay 
from high-spin states induced by heavy-ion compound reactions is shown in 
fig. 19. 

The 2-D and 3-D fluctuation analysis method developed by this group 
[55] has been applied to the spectra from a triple coincidence experiment 
made by the Manchester-Daresbury-Copenhagen-Bonn collaboration at the 
Daresbury Tandem Accelerator. The reaction consisted in bombarding 124 Sn 
by a beam of 4S Ca of 215 MeV, to form the residual nuclei of 168 Yb at the 
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highest spin I max = 60h. From the fluctuation analysis done on 168 Yb they 
concluded that the rotational correlations originating from well defined ro- 
tational band structures only exist up to 1 MeV above the yrast line. Over 
this limit the main decay routes spread out into many branches. This may 
be accounted for by the rotational damping phenomenon, as Guhr and Wei- 
denmuller [56] had previously suggested in 1989 from a theoretical calcula- 
tion with a GOE hamiltonian weakly perturbed by a residual interaction. In 
the conclusions of this work, the authors proposed the measurement of the 
spreading width of E2-strength for a high-spin state well above the yrast 
line, to determine the possible onset of chaos in "warm" deformed nuclei. 
The fluctuation analysis mentioned above has in fact taken up this sugges- 
tion. 

The first conclusions extracted from this new analysis of 7-ray data sug- 
gest that there is much interest to be had from the study of states near 
and above the yrast line, "warm" rotational nuclei and the decay of unsta- 
ble quantum states. New theoretical attention has been devoted to these 
phenomena. 

6.3 Rotational damping motion 

Due to the recent experiments in 7-ray detection, a renewed interest has 
developed in the the study of the decay of chaotic states and the possibility of 
clarifying the mechanisms by which an excited nucleus (6-8 MeV over yrast 
line) undergoes a transition from a regular to a chaotic behaviour. 

In this sense, the work of Matsuo, Dossing, Herskind and Frauendorf [57] 
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is very interesting. They investigate the properties of energy levels and rota- 
tional transitions as a function of the excitation energy, in "warm" deformed 
nuclei, and show that, at an excitation energy of the order of E x ~ 8 MeV, 
there is considerable fragmentation of the rotational E2 strength distribution 
when it is represented against the 7 energy (see fig. 20). This fragmentation 
with increasing energy is due to the mixing of the rotational bands caused 
by the Surface Delta Interaction (SDI): 

1/(1,2) = AnVoSin - R )5(r 2 - R ) • E^*>i) • ^(r 2 ), (6.6) 

which contains all the possible multipolarities. Actually, the high multipole 
components of this SDI are shown to be responsible for the mixing of rota- 
tional bands which is reflected in fluctuations typical of chaotic behaviour. 
The strength fluctuations of the E2 transitions, for instance, obey a Porter- 
Thomas distribution above a certain excitation energy, assumed to be the 
threshold for the realization of quantum chaos in the system (fig. 21). 

In particular, the same authors have also studied high spin levels in nuclei 
with a cranked shell model extended to include residual two-body interactions 
[58] . Here again the residual interaction induces the transition from a regular 
to a chaotic regime, since when the rotational bands do not interact, 7-ray 
energies behave like random variables, i.e., they obey a Poisson distribution 
typical of a regular regime. On the contrary, when the residual interaction is 
added, at an excitation energy over 600 KeV, a gradual rotational damping 
is established, and, at 1.8 MeV above the yrast line, the complete damping 
is observed and typical GOE fluctuations of the energy levels and transition 
strengths are produced (fig. 22). 
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A similar approach is used by Aberg [59] to study the rotational damping 
of rapidly rotating nuclei. As in the preceeding approach, a cranked Nilsson 
model with a schematic residual interaction is used to study the distribution 
of the E2-strength, and 7 — 7 correlations. In particular the normal deformed 
(ND) l68 Yb and the superdeformed (SD) 152 Dy are studied, and the different 
behaviours of these nuclei with increasing energy are compared. The main 
difference between the approach of Aberg and that of the authors of ref. 
[58] is that in the former the energies of the np-nh states are given in the 
laboratory frame as functions of the angular momentum I. So, diagonalization 
is performed for each fixed I value. This is reflected in a change of the 
individual energies, but the conclusions regarding the different statistical 
properties are the same. The Aberg approach is perhaps more useful if a 
direct comparison with E 1 —E 1 experimental spectra is to be carried out. The 
"classical" eigenvalue statistics are studied (i.e. -P(s), A 3 (L)) with special 
regard to the A 3 (L) since it gives a measure of long-range correlations. In this 
model the transition from order to chaos is brought about by changing the 
strength of the two-body force (A) which is taken as a parameter, V 2 = ±A, 
where the sign is chosen randomly to avoid coherent effects. Following the 
procedure of Brody [60], the A 3 (L) spectrum is fitted with a single parameter, 
q, to study the mixing between the Poisson and GOE distribution: 

A 3 (L, q) = Af [(1 - q)L] + A° OE (qL). (6.7) 

The dynamical properties are studied in this model and again a rotational 
damping is manifested in the fragmentation of E2 strength in many daughter 
states. The average standard deviation of the E2-strength function saturates 
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for sufficiently large values of the strength of the two-body residual interac- 
tion, following the same trend as q parameter, which gives the mixing of a 
Poisson and GOE behaviour of spectrum properties (see figs. 23-24). 

In addition to the previously mentioned rotational damping, Aberg also 
studies motional narrowing, i.e., the decrease of the width of nuclear mag- 
netic resonances when the temperature increases, a phenomenon which may 
be studied in a simple two-band model. The narrowing of the strength func- 
tion is accomplished by a change from a Gaussian shaped strength to a 
Breit-Wigner shape. This phenomenon is understood in terms of time scales: 
when the time scale of the fluctuations in available rotational frequencies is 
much less than the time the wave function takes to spread out over the basis 
states, the intensity spectrum is Gaussian, while in the opposite case the 
corresponding spectrum is of the Breit-Wigner type. 

Since, due to the motional narrowing the E2 strength may become very 
narrow at the end, a small number of relatively strong E2 transitions may be 
observed, and, consequently, a " suppression of chaos" . In fact the transition 
entropy defined as [44]: 

H T " = -^[M ae ,] 2 -ln([M ao ,] 2 ), (6.8) 

a' 

where: 

M aa ,(I)=<a',I-2\M(E2)\a,I> (6.9) 

decreases when motional narrowing sets in. 

The 7 — 7 correlations between consecutive 7-rays are also analyzed to 
study the fragmentation of E2 strength. Correlations are revealed at small 
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excitation energies or when a small A is used, while at high excitation energies 
these correlations disappear. 

6.4 The study of unstable systems using a complex effective 
hamiltonian 

A different and more general approach to the study of unstable quan- 
tum systems is that of Sokolov and Zelevinsky [61] and more recently that 
of Mizutori and Zelevinsky [62], in which the statistical theory of spectra, 
formulated in terms of random matrix theory, is generalized to treat unsta- 
ble states, i.e. those coupled to open channels. In this way the influence of 
the coupling with the continuum on the properties of internal states can be 
better understood. 

In general, as we have seen up to now, the level statistics are treated as if 
the states were stationary but, in reality, all excited states are resonances em- 
bedded in the continuum. So if one intends to study the most general states 
of a system (including the unstable ones) a new approach to the problem has 
to be used. 

When the widths of the states are small as compared with level spac- 
ings, the approximation of discrete levels is reasonable for long-lived states. 
However, when the widths increase and the levels overlap, the effects of the 
finite lifetime of those become crucial and the application of a random matrix 
ensemble of hermitian hamiltonians is not enough appropriate. 

Since the standard gaussian ensembles are applicable to discrete station- 
ary states, while an excited state decays via open channels, the use of a 
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non-hermitian effective hamiltonian H is imposed to take into account the 
width of the states due to their finite life [61,62]. 

Within this model the reaction amplitudes are represented, using the 
general theory of resonant nuclear reactions, as sums of pole terms in the 
complex plane. Such poles, e n = E n — {l/2)iT n are the eigenvalues of the 
non-hermitian effective hamiltonian H = H — (i/2)W where H is the hermi- 
tian part belonging to the GOE. The amplitudes of the antihermitian part, 
W, have a separable structure W nm = Yi c A c n A c m due to the unitarity of the 
scattering matrix, where A c n are the amplitudes for the decay of intrinsic 
states of the system, \n > (n=l,...N) into different channels, c (c=l,....k) 
and that must be real according to time reversal invariance of the full hamil- 
tonian, H. Those eigenvalues correspond to the intermediate unstable states 
of energies E n and widths Y n . 

Decaying systems are hence described by ensembles of random non-hermitian 
hamiltonians represented by N-dimensional matrices where the hermitian 
part of the effective hamiltonian, H, belongs to the GOE. The decay ampli- 
tudes A^, on the other hand, are assumed to be Gaussian random variables 
completely uncorrelated with the matrix elements of H and among each 
other. 

Within this model the distribution function of complex energies, the level 
spacing distribution and the width distribution are studied in ref. [61,62] for 
weak and strong coupling to the continuum, as well as for the transitional 
region. 

The authors of ref. [62] discussed in great detail the single channel case: 
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(6.10) 



The strength of coupling to the continuum is given by the overlap parameter 
X = rj/a = 2<T> / D where rj — < r? >; D is the mean level spacing. 

When the overlap parameter is small (x << 1), that is, when coupling 
with the continuum is weak, the hermitian part H dominates, while the 
antihermitian part, W, prevails in the opposite limit (x >> 1). 

Typical signatures of the transition from weak to strong coupling are ob- 
served: first, concerning the level spacing distribution, which gives the sim- 
plest characterization of spectral correlations, it is observed that for a weak 
coupling (where the hermitian part, H, of the effective hamiltonian domi- 
nates) the distribution of levels corresponds to a Wigner distribution, and 
deviations from it are small. As the coupling with the continuum increases, 
the Wigner function ceases to be a good approximation to the GOE and a 
new feature is observed: the level repulsion at short distances disappears. 

The authors of ref. [62] have fitted the P{s) distribution by a simple 
one-parameter superposition of the normalized Wigner, P w , and Gaussian, 
Pg distributions, that is: 



where the parameter A is determined, for a fixed a, by the mean level spacing: 



In the fig. 25 the trend of the mixing parameter a for the P(s) as a func- 
tion of the coupling constant x is plotted, showing a pronounced minimum 



la 2 s 2 

P(s) = aP w (s) + (1 - a)P G (s) = -[-s + (l- a )J-] exp (-_), (6.11) 




(6.12) 
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(a = 0.74) at the transition point x — 1- Both small and large values of x 
correspond to the nearest level statistics close to a Wigner distribution P w . 
It is well known that the GOE implies rigid spectra with small level spacing 
fluctuations. So, we can conclude that coupling to the decay channels softens 
the spacing distribution so that the spacing fluctuations increase gradually 
in the transition region. 

As regards the distribution of complex energies, it is has been observed 
that, for a random matrix ensemble, the widths distribution corresponds to 
a Porter- Thomas one, i.e.: 

P PT (r) = . L— exp (-^r). (6.13) 

y 2,71 < r > /r) 2 

This, however, is no longer correct when the level spacings become very small 
compared to the widths, i.e., when the overlap between neighbour states is 
considerable. This is a consequence of the coupling to the continuum, since 
for weak coupling the typical widths are small compared to spacings. When, 
however, the coupling to the continuum is important, the widths become 
larger than the mean level spacing. This behaviour had already been observed 
for complex random matrix ensembles [63] and in chaotic dissipative systems 
[64]. 

The most interesting statistic for highlighting the phase transition is the 
width distribution. In fig. 26 the widths, plotted as a function of the coupling 
parameter x, show a collectivization. The mean fraction < Ti/^Fj > of 
the total width is accumulated by the broadest state IY 

In conclusion, we can say that within these works a standard statistical 
spectroscopy of discrete levels and unstable states, as well as for the transi- 
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tional region, has been performed using ensembles of random non-hermitian 
hamiltonians represented by N-dimensional matrices. They show [61,62] that 
the instability of states changes the statistics remarkably, removing level re- 
pulsion when distances are smaller than widths. When the coupling to the 
continuum increases, that is, when the matrix elements of the antihermitian 
part of the effective hamiltonian become comparable with the spacings of 
eigenvalues of the hermitian part, a transition to a new regime occurs. 

Related to the study of unstable systems and in the same spirit as the 
preceeding work, a paper by Haake et al. [65] is noteworthy. The authors 
study the level density of different classes of random non-hermitian matrices, 
H = H + iT, where the damping, T, is chosen quadratic in Gaussian random 
numbers, to describe the decay of resonances through various channels. When 
the notion of level spacing is extended to the Euclidian distance between com- 
plex eigenvalues of non-hermitian operators, two different behaviours of P(s) 
statistic (as in the ordinary case of real eigenvalues with the Wigner and Pois- 
son statistics) allow the distinction between regular and chaotic dynamics: 
cubic repulsion tends to be typical under conditions of global classical chaos 
while linear repulsion signals classical integrability. This classification also 
arises for dissipative systems [66]. 

6.5 The IBM model in the study of unstable systems 

The profound understanding of the statistics of low-lying levels of nuclei 
and the underlying chaotic dynamics requires realistic theoretical models of 
the nucleus. In general, the models used have only two degrees of freedom. 
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However, since the quadrupole deformation plays an important role in col- 
lective nuclear dynamics, a realistic model requires at least five degrees of 
freedom. 

The work of Whelan and Alhassid [67] have been along these lines. Us- 
ing the Interacting Boson Model (IBM) [3], whose classical limit is obtained 
with coherent states and which provides a good description of low-lying en- 
ergy levels and EM transitions of heavy nuclei, they studied, classically and 
quantum mechanically, the chaotic properties of low-lying collective states of 
atomic nuclei. The six degrees of freedom come from one monopole s-boson 
(0 + ) and five (/x = —2,..., 2) bosons (2 + ) with which a U(6) dynamical 
algebra is constructed. The most general hamiltonian is then constructed 
with all one and two-body scalars that conserve the total number of bosons, 
N = ss + + Z li d+d li . 

The most useful parameterization of the IBM hamiltonian is given by: 

H = E + c n d + c 2 Q x -Q x + c l L 2 , (6.14) 

where n d = d + ■ d is the number of d-bosons, L is the angular momentum 
and Q x is the quadrupole operator: 

= (d + x s + s + x d) {2) + X {d + x d) {2 \ (6.15) 

that depends on a parameter x, and where d^ = (-)^d^ so that d M transforms 
under rotations like dt- 

Its classical limit is obtained for the number of bosons going to infinity, 
since 1/N plays the role of K. That classical limit is: 

h = e + c[r]n d - (1 - r])q x ■ q x ] + CiL 2 , (6.16) 
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where: 
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rj 



(6.17) 
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with < rj < 1. In an algebraic model like this, a dynamical symmetry 
occurs when the hamiltonian can be written as a function of the Casimir 
invariants C of a chain of subalgebras of the original algebra. So, the authors 
of ref. [67] analyze the rotational nuclei (SU(3) symmetry), vibrational nuclei 
(U(5) symmetry) and 7-unstable nuclei (0(6) symmetry): 

U(6)DU(5)DO(5)DO(3) (vibrational nuclei) 

U(6)dSU(3)dO(3) (rotational nuclei) 

U(6)dO(6)dO(5)dO(3) (7-unstable nuclei). 

The authors study the character of the classical dynamics of a nucleus 
described by the classical limit of the quantum IBM hamiltonian. As a first 
result, they found that the quantal fluctuations, which are well correlated 
with the classical results, are independent of the number of bosons. The 
transitions between the different dynamical symmetries of the model (i.e., 
between the different dynamics of nuclei) with the variation of the parameter 
X, of both classical and quantal hamiltonian, are studied. 

The transition between deformed rotational nuclei (SU(3)) and spherical 
vibrational nuclei (U(5)) is obtained for \ — ~ 1/2a/7 and < r\ < 1. With 
these parameters, classical and quantal measures of chaos are studied: the 
average maximal Lyapunov exponent A and the fraction of chaotic trajecto- 
ries a, on the one hand, and the parameter uj of the Brody distribution of 
level spacing [60] on the other: 



Pw(s) = As u exp (—as' 



)■ 



(6.18) 
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The Brody distribution, as the A|(L) mentioned above, interpolates between 
a Poisson distribution (cu = 0) for a regular system and a Wigner one (cu — 1) 
for a chaotic system, a and A are chosen such that P w (s) is normalized to 1 
and < S >= 1. A second quantal measure of chaos is v>, which characterizes 
the B(E2 :/—>/) distribution for the levels with I = 2 + : 

p " fo) = ftfyfe) 5 ^ exp( -2777»- (619) 

where: 

2 /=|</|TK>| 2 , (6.20) 

|i > and \f > being, respectively, the initial and final states and T the 
transition operator. P(y)dy is the probability of having intensity P(y) in 
the interval dy around y. P(y) reduces to a Porter-Thomas distribution for 
v — 1 (chaotic limit) and, as the system becomes more regular, v decreases 
towards small positive values. 

The maximum of chaos is obtained for rj = 0.5 — 0.7. The quantal re- 
sults show a strong correlation with the classical dynamics, since uj and v 
are largest around r\ = 0.5 — 0.7. Secondly, the transition between rota- 
tional nuclei (SU(3)) and 7-unstable nuclei (0(6)) is obtained for rj = and 
— l/2y/7 < x < 0, chaos being settled for intermediate values of x- Finally, 
the transition between 7-unstable (0(6)) and vibrational nuclei (U(5)), for 
X = and < 77 < 1, is always completely regular. 

The most important result is, on the one hand, the discovery of a new 
nearly regular region that lies between a rotational (SU(3)) and vibrational 
(U(5)) regime of the nucleus, and that is not related to any of the known 
dynamical symmetries of the model. Since the fraction of chaotic trajectories 
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(er) is < 0.3, that region is not completely regular. The signatures of this 
nearly regular region are identified by the rather sharp minimum in all mea- 
sures of chaos, as seen in fig. 27. The authors suggest that this new nearly 
regular region may be connected with a previously unknown approximate 
symmetry of the model. On the other hand, a spin dependence of the degree 
of chaoticity is revealed. When both classical (A,<r) and quantal (uj,v,q) mea- 
sures of chaos are plotted versus spin, an interesting dependence is found: in 
all measures (except A) a weak dependence on the spin, at low and medium 
spins, is found. However, at high spins (I > 20h) there is a rapid decrease of 
chaoticity and the motion becomes regular. 

In conclusion, strong correlations are observed between the onset of clas- 
sical and "quantum" chaos. Although the IBM model can be useful for the 
study of the degree of chaoticity of low-lying collective states of nuclei, at 
higher spin and/or energy, bosons may break into quasi-particles and it is also 
important to take into account the additional fermionic degrees of freedom, 
to obtain a realistic description [4]. 

In this section we have presented some of the most interesting recent 
works, which have attempted to explain the mechanisms of the decay of 
quantum states and its connection with the onset of chaos. 

The theoretical efforts are also supported by the new detector genera- 
tion, which, with their higher resolution, can provide a new and profound 
knowledge in the field of nuclear spectroscopy. 
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7. Conclusions and open problems 

To conclude this presentation of some regular and chaotic aspects in nu- 
clear dynamics, it is important to remember two points. The first is the good 
agreement between the theoretical previsions concerning the Poisson-^GOE 
transition and the experimental data; the second is the possibility of studying 
the order to chaos transition using simple models. 

As pointed out by many authors (see for instance ref. [1]), the static 
nuclear mean field is too regular (as witnessed by the existence of shell struc- 
ture) to be responsible of chaotic behaviour, and chaos must be caused by 
residual interaction, as shown in the models discussed above. Because of the 
exclusion principle, the role of the latter increases with the excitation energy. 
This makes one wonder whether nuclear motion might be more regular in the 
ground-state region than, say, at neutron threshold. Indeed, at several MeV 
excitation energy, collective states and giant resonances acquire substantial 
spreading widths [1]. They lose their individuality as excitation states, and 
show up only as bumps in strength functions. The same kind of statement 
applies also to rotational bands 1 or 2 MeV above the yrast line. Much work 
remains to be done to elucidate the role of chaotic motion and its interplay 
with collectivity. 

There is also the problem of showing which properties of the nucleon- 
nucleon interaction justify the use of random matrix models. These, accord- 
ing to ref. [1], are related to nonintegrability and chaos, but the matter 
requires further research to obtain a deeper understanding of the problem. 

We have shown that chaos in nuclei is a rich field; it permeates many 
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aspects of nuclear structure and reactions, and gives rise to new analyses of 
the experimental data. 
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Figure Captions 

Figure 1: The total neutron cross-section on 232 Th vs neutron energy 
(adapted from [7]). 

Figure 2: The NDE experimental results compared with the theoretical 
predictions (adapted from [19]). 

Figure 3: P(s) for "cold" deformed rare-earth nuclei (adapted form 
[20]). 

Figure 4: Spectral statistics for nuclei with 24 < A < 244 and excitation 
energy of few MeV; (a) 2 + and 4 + states, even-even nuclei; (b) all other 
states, even-even nuclei; (c) states with non-natural parity, odd-odd nuclei; 
(d) states with natural parity, odd-odd nuclei (adapted from [21]). 

Figure 5: A diagramatic representation of the ranges of excitation energy 
E x and angular momentum / associated with each set of data previously 
discussed (adapted form [20]). 

Figure 6: Closed (a) and not closed (b) orbits on the torus. 

Figure 7: A Poincare section. 

Figure 8: Poincare sections: (a) regular closed orbit, (b) regular not 
closed orbit, (c) chaotic orbit. 

Figure 9: P(s) and A 3 (L) for billiards. Regular billiard: circular bound- 
ary (left), chaotic billiard: stadium boundary corresponding to eigenfunctions 
with odd-odd symmetry (right) (adapted from [29]). 

Figure 10: E — T plots for different values of x and for different families: 
A initial conditions along the axis q±, □ initial conditions along the axis q2, 
O initial conditions near the minima of the static potential: (a) x — 0-5, (b) 
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x = 2,(c)x= 100 (adapted from [43]). 

Figure 11: Ratio between the number of stable periodic orbits and the 
number of total periodic orbits vs x (adapted from [45]). 

Figure 12: Universal and long-range behaviour of A 3 (L) vs L. The 
continuous curve is the GOE behaviour (adapted from [43]). 

Figure 13: P(s) vs s for the eo class with M = 102; (a) x = 0.75, (b) 
X = 2, (c) x — 3. The continuous curve is the Brody distribution (6.12) with 
u = l- AN st /AN tot (adapted from [45]). 

Figure 14: Density of quasi-crossing vs x for all the classes with M = 
102 (adapted from [45]). 

Figure 15: A 2 (E) vs energy E for different values of x f° r the eo class 
with M = 102 (adapted from [45]). 

Figure 16: Momenta I m vs energy E for the eo class with M = 102: (a) 
chaotic region, (b) quasi-integrable region (adapted from [51]). 

Figure 17: Information measure S vs energy E for the eo class with 
M = 102: (a) chaotic region, (b) quasi-integrable region (adapted from 
[51]). 

Figure 18: Transition probability (ee) — > (oo) with M = 50: (a) regular 
region, (b) chaotic region (adapted from [53]). 

Figure 19: A schematic illustration of the average flow of 7-ray decay 
from high-spin states induced by heavy-ion compound reactions. The insert 
illustrates the spreading of the E2 transitions (adapted from [55]). 

Figure 20: The rotational E2 strength distribution vs the energy E 1 
for four initial levels with different excitation energies above the yrast line 
(adapted from [57]). 
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Figure 21: The probability distribution of the rotational E2 strengths 
Sij for 7-ray energies in the interval 0.90 < E 1 < 1.05 MeV. The dashed curve 
represents the Porter-Thomas distribution. The excitation energy above the 
yrast line is indicated for each bin (adapted from [57]). 

Figure 22: The rotational-strength distributions obtained for transitions 
from 50 levels in each bin. The excitation energy above the yrast line is 
indicated for each bin. The smoothed distribution function is drawn with a 
solid curve. The transitions shown are for (— , 1) in 168 Yb at Tiw = 0.5 MeV 
(adapted from [58]). 

Figure 23: Mixing parameter q vs the strength of the two-body inter- 
action A, at two excitation energy intervals, 2.0 — 2.5 MeV (dashed line) 
and 3.0 — 3.5 MeV (solid line) in superdeformed 152 Dy, and 0.5 — 1.0 MeV 
(dashed line) and 1.5 — 2.0 MeV (solid line) in normal-deformed 168 Yb. In 
all cases I w = 50 + (adapted from [59]). 

Figure 24: Average standard deviation of the E2-strength function vs 
A for 50 + — > 48 + transitions emerging from the excitation energy intervals 
2.0 — 2.5 MeV (dashed line) and 3.0 — 3.5 MeV (solid line) in superdeformed 
152 Dy, and 0.5 - 1.0 MeV (dashed line) and 1.5 - 2.0 MeV (solid line) in 
normal-deformed 168 Yb (adapted from [59]). 

Figure 25: Mixing parameter a vs the continuum coupling constant \ 
(adapted from [62]). 

Figure 26: Mean ratio of the width of the broadest state to the total 
summed width of all complex eigenvalues vs the continuum coupling constant 
X- Error bars correspond to the spread over 50 matrices 160x160 (adapted 
from [62]). 
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Figure 27: Classical and quantal measures of chaos vs v (v = 0: SU(3), 
v — 1: U(5)). Left: A and a for angular momentum per boson / = 0.1. 
Right: oo and v for the levels with I = 2 + (adapted from [67]). 
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